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We study the location of the partition function zeros in the complex )3 plane (Fisher's Zeros) for 
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L = 4 and 6. We compare results obtained using various polynomial approximations of the loga- 
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at different j3 provides the more reliable estimates. 
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1. Introduction 

Locating the zeros of the partition function of lattice gauge theories in the complex j(3 plane and 
their volume dependence is important to understand the large order behavior of the weak coupling 
expansion [|l|, ^ |3[ ^ at zero temperature and the nature of the finite temperature transition |6|]. 
These zeros are called Fisher's zeros [^] and should not be confused with Lee- Yang zeros which 
are zeros in the complex fugacity plane or the complex exp(— 2j8//) plane [Q] . 

In the following, we discuss the Fisher's zeros of a pure gauge SU (2) theory with a partition 
function 

Z = Y{f dUit-^^ , (1.1) 

where S is the Wilson action 

S = Y.i^-{imReTr{Up)), (1.2) 
p 

and j8 = 2N/g^. 

Our expectations is that at zero temperature, there is no singularity on the real axis of the 
complex j8 plane and as the volume increases, the zeros stay away from the real axis. On the other 
hand at non-zero temperature, we expect that as the volume increases, the zeros pinch the real axis 
as for the 2D Ising model |^]. 

The zeros of the partition function can be calculated using the reweig hting method [|, |]. 

Z(i3o + Aj8) =Z(A)) < exp(-Ai35) >ft . (1.3) 

It is convenient to subtract < S >^ from S in the exponential because it removes fast oscillations 
without changing the complex zeros. 

Z(j8) is the Laplace transform of the density of states n{S): 

Z{I5) = / dSn{S)exp{-l5S) (1.4) 

JO 

One can show [ jl0| ] that for SU (2) Smax = 2^ and that for lattices with an even number of sites 
in each direction, n{S) = n{2.yVp — S), where ,yVp is the number of plaquettes. For a D dimensional 
cubic lattice with periodic boundary conditions, jVp = LPD{D— 1)/2. When n{S) is known, it is 
possible to calculate the partition function for any complex value of j8. 

In the following, we consider symmetric 4^ and 6"* lattices. For values of j3 near 2, the dis- 
tribution of S is nearly Gaussian and the location of the peak scales with the number of sites. The 
departure from a Gaussian distribution is hardly visible on a histogram. However, as shown in Fig. 
[T], the residuals show a coherent behavior on a 4^ lattice. As the volume increases, the non-Gaussian 
features are scaled down and for j3 = 2.18 it seems that the signal is lost in the statistical noise. In 
this figure, Ni is the number of data points in the /-th bin and Pj the corresponding probability for 
a Gaussian distribution with the estimated mean and variance. As in the Gaussian approximation 
there are no complex zeros. It is crucial to resolve the departure from this approximation. We 
now discuss two methods, one based on the estimation of the moments and the other on numerical 
calculation of the the density of states. 
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Figure 1: The residuals {Ni-NPi)/{NPi)^l^ discussed in the text for a distribution of 500,000 values of S 
in an histogram with 100 bins for a SU{2) pure gauge theory on a 4"^ lattice at j3 = 2.18. 

2. The Moments Method 

In this section, we consider the following corrections |^ to the Gaussian approximation: 

P{S) oc exp(-Ai5 - X2S^ - X3S^ - X^S"^) (2.1) 

The four unknown parameters can be determined from the first four moments using Newton's 
method. The moments are defined as 

mi = <S> /J^p 
m2 = <{S-<S >)2 > /^p 
= <{S-<S >)^ > /^p 

1714 = (< {s- < 5 >)^ > -3 < {s- < s >)2 >^)/-yyp 

(2.2) 

As < 5 > scales like and the individual terms of the n-th moment like each subtraction 
imphes a loss of significant digits which increases with the volume. As shown in Fig. ^ the third 
and fourth moments have large errors even on a 4^* lattice. Once we obtain Z(j3), we can calculate 
the zeros of real and imaginary parts of it separately. The cross points are the Fisher's zeros. The 
result for jSq = 2.18 on a 4^ and 6^ lattice is shown in Fig. ^ The errors of the moments affect 
P{S) and the location of the zeros. A change of the fourth moment within the error bars produces 
changes in the zeros illustrated in Fig. ^. This change gives an idea of the errors associated with 
the method. 



3. Density of state Method 

The probability distribution of the plaquette can be written as 

Pp{S)=n{S)exp{-l5S) . (3.1) 
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Figure 2: and niA, for SU(2) as a function of j3 on 4 and 6 lattices 



To find the j8 independent density of state using Monte Carlo data, we needed to patch the data 
from different jS together. First the j3 dependence was removed by multiplying by e^^ . Using only 
the bins with statistics higher than half the maximum, we overlay the data from each set on top 
of one another to make a smooth curve (we took the log of the values in the bins and adjusted 



the offset with a one-parameter fit). This procedure can be found with more detail in Using 
numerical interpolation for f{S/^p) = ln{n{S / J\^p)) / -yVp, it is possible to calculate the zeros using 
numerical integration. The results are shown in Fig. ^. 

As the changes were more important than expected (compared to Fig. Hb, we estimated the 
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Figure 3: Zeros of the real (crosses) and imaginary (circles) using MC on a 4^ lattice, for5t/(2) at j3o = 2.18. 
The smaller dots are the values for the real (green) and imaginary (blue) parts obtained from the 4 parameter 
model. The MC exclusion region boundary for d — 0.20 (defined in [|[]) is represented by boxes (red). 



errors by using different distributions of S. Three different methods are used to get the partition 
function, with which < cos(/mj8(5'— < S >)) > is calculated and compared as a function of the 
imaginary part of j8 at fixed real part 2.18. We estimated the error by generating multiple data. For 
the first, we generated 50 bootstraped sets of configurations and computed < cos{Imp{S— < S > 
)) > directly by MC average. For the second, we get the partition functions via the density of states 
which are obtained using interpolation out of 50 patchings. For the last, we fit the 50 patchings 
using Chebyshev Polynomials instead of interpolation. The Chebyshev fitting seems to have much 
higher accuracy and stability than the other two methods and will be used for further investigations. 
Calculations of the zeros for Imp < 0.11 with this method appear to be consistent with Fig. ^. 
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Figure 4: Changes of the zeros when the fourth moment is shifted in the positive direction until the error 
bar is reached for j3o = 2.18 on a 4^ lattice. 
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Figure 5: Same quantities as in Fig. 4 on a 4^ lattice but with an interpolated version of /. 
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Figure 6: < cos(/mj3(5— < S >)) > as a function of the imaginary part of j3 at fixed real part 2.18 with 
three methods: interpolation, Chebyshev fitting and bootstraped Monte Carlo, are on 4"^ (left) and 6'^ (right) 
lattices . 
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4. Conclusions 

In conclusion, we have compared the moments methods and methods based on the density of 
states with the simple MC reweighting procedure to calculate the zeros of the partition function. 
The method where the density of states is approximated by Chebyshev polynomials seems the most 
reliable and will be used in future investigations. 
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